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Abstract 

We present an analysis of the bp production from double parton (DP) sact- 
tering and single parton (SP) scattering in the large px region via color-octet 
gluon fragmentation. We find that at the Tevatron the DP pp production 
is at the edge of the detectability at present, and at the LHC the DP cross 
section will dominate over the SP cross section in the lower pr(min) region 
(he., p T (min) < 7 GeV). We also conclude that the color-octet mechanism 
is of crucial importance to the double J/ij) production at high energy hadron 
colliders. 

PACS number(s): 13.85.Ni, 14.40.Gx 
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The multiple parton (MP) interactions, in which two or more pairs of partons scatter off 
each other, might become important parts of the structure of typical collisions at very high 
energy hadron colliders, mostly because of the rapidly increase of the cross section for pro¬ 
duction of jets with transverse momentum px > px(nnn) ps 2 GeV. Also, the Double Parton 
scattering (DP) processes, where two parton-parton hard scattering takes place within one 
pp (or pp) collision, have always been a topic to study possible parton-parton correlations 
in the proton |l[. In experiment, to search for the DP signals, one must subtract the back¬ 
ground signals from the single parton (SP) scattering. Any possible process to study the DP 
signals must have a large enough ratio of signal to background events. This is important to 
theoretical investigations. The mostly considered DP process is the production of four liigli- 
Pt jets via double parton scattering within a single hadronic collision J^] (4 —> 4 reactions). 
In such process, two pairs of jets are produced in the final states and each pair has equal and 
opposite transverse momentum. Various hadron collider experiments have searched for this 
signature |3|-|5|. However, the four jets production from the SP processes (2 —> 4 reactions) 
has a large background. For example, the CDF collaboration at the Fermilab Tevatron 
found evidence that 4 —» 4 reactions contribute only about 5% to the production of four 
jets with px > 25 GeV [(|. Another sort of the DP signals is the double Drell-Yan processes 
[(|, in which double lepton-pair production may provide a cleaner signature than that of 
four-jets process. However, such processes which are due to quark and anti-quark annihila¬ 
tions are much less important compared to gluon-induced processes in higher energy hadron 
collisions. More recently, the production of one-photon + three-jets as signals for the DP 
interactions has been studied [0, and a strong signal is observed by the CDF collaboration 
at the Fermilab Tevatron ||. 

In this note, we discuss the probability of double J/pj production at large px as signals 
for DP scattering. J/p) production at hadron colliders is of special significance because it has 
extremely clean signature through its leptonic decay modes. In recent years, J/pJ production 
at large px is well studied both in theoretical and experimental sectors. J/ip production 
at large px has two main sources. One is from b decays, and the other is from so-called 
“prompt” production. The contributions from b decays can be isolated by using a secondary 
vertex detector, then the prompt J/pJ production can be studied in experiment [tj. At high 
energy hadron colliders, unlike the Drell-Yan production mostly coming from the quark and 
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antiquark annihilations, prompt J /ip production dominantly comes from gluon-gluon fusion 
process (gg —*• ipX). This gluon-induced process would provide much larger cross section 
for the production, and then can be used to study other interesting physics, such as the DP 
scattering process. By triggering four-// final states to measure the double J/ip production 
at large px, we can search for the DP signals in high energy hadronic collisions. 

In PU[ , the authors have calculated ipip production by considering the single J/ip produc¬ 
tion at the leading-order within the conventional Color-Singlet Model. However, as pointed 
out by Braaten and Yuan (ll|], the fragmentation contributions will dominate over those from 
leading-order processes at sufficiently large px, although the fragmentation processes are of 
higher order in strong coupling constant a s . Explicit calculations of the contributions to ip 
production at the Tevatron from fragmentation of gluons and charm quarks revealed that 
fragmentation dominates over the leading-order gluon-gluon fusion mechanism for px greater 
than 6 GeV 0. More recently, the CDF collaboration have reported their measurement of 
charmonia production at large px- They found a large excess of direct production (excluding 
the contribution from b decays and the feeddown from y c ) both for J/ip and ip' 0 [|T3|j . 
The experimental measurement is a factor of 30 ~ 50 larger than the theoretical prediction 
of the Color-Singlet Model. Motivated by this “surplus” problem, a new mechanism for 
heavy quarkonium production at large px, named as Color-Octet gluon fragmentation has 
been proposed [|T|], which is based on the factorization formalism of non-relativistic quan¬ 
tum chromodynamics (NRQCD) |T5|. In this approach, the production process is factorized 
into short and long distance parts, while the latter is associated with the nonperturbative 
matrix elements of four-fermion operators. This factorization formalism provides a new pro¬ 
duction mechanism, the color-octet mechanism, in which the heavy quark and antiquark 
pair is produced at short distance in a color-octet configuration and subsequently evolves 
nonperturbatively into physical quarkonium state. In the past few years, applications of the 
NRQCD factorization formalism to J/ip{ip') production at various experimental facilities 
have been studied H- 


Single J/ip production at large px may dominantly come from gluon fragmentation con¬ 
tributions. According to NRQCD factorization formalism, the gluon fragmentation to J/ip 
production can be factorized as, 




( 1 ) 
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where z is the longitudinal momentum fraction carried by the produced J/i[} in gluon frag¬ 
mentation, n = 2 m c is the fragmentation scale. d g ^ n represent the short-distance coefficients 
associated with the perturbative subprocesses in which a cc pair is produced in a configu¬ 
ration denoted by n (angular momentum 2S+1 Lj and color index 1 or 8). ( O j J are the 
long distance nonperturbative matrix elements demonstrating the probability of a cc pair 
evolving into the physical state J/ij). The short-distance coefficients d g ^ n can be obtained 
from perturbative calculations in powers of coupling constant a s . consist of two 

kinds of matrix elements, he., the color-singlet and color-octet matrix elements (according 
to the color index 1 or 8). The color-singlet matrix elements may be related to the quarko- 
nium radial wave function or its derivatives at the origin, and may be calculated by potential 
models or estimated by leptonic decay widths of quarkonium states. Whereas the color-octet 
matrix elements can only be determined by fitting the theoretical prediction of quarkonium 
production rates to the experimental data. The relative sizes of these matrix elements can 
be determined by their scale properties with v 2 according to the NRQCD velocity scaling 
rules, where v is the typical relative velocity of the heavy quark in the bound state. The 
fragmentation function in Eq.(lj) is a double expansion in a s and v. 

For J/iJj production in gluon fragmentation, the color-octet matrix element (O'l^ ( 3 S\)) 
is smaller than the color-singlet matrix element ( 3 <S'i)) by a factor of order v 4 according 

to the NRQCD velocity scaling rules. However, the short-distance coefficient for the color- 
octet term in Eq.(JiT) is larger than that for the color-singlet term by a factor of order 1/ol 2 s . 
Numerical results show that color-octet contributions are 50 times larger than color-singlet 
contributions ]h||. In the following calculations, we neglect the color-singlet term in gluon 
fragmentation in Eq. ([!]), and only consider the color-octet gluon fragmentation. The leading- 
order color-octet gluon fragmentation to J/tp production gives JT4J] 

( 2 ) 
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In our calculations, the effects of the evolution of gluon fragmentation function with scale 


H 1 are neglected, which may introduce some error. However, as argued in [|lj], including 
evolution would not necessarily be an improvement, since naive Altarelli-Parisi equations do 
not respect the phase-space constraint D g ->j/^(z, /i 2 ) = 0 for z < m 2 ,^/p 2 [|T^ . 

We plot in Fig.l the typical Feynman diagrams for double J/ijj production at large 
Pt via color-octet gluon fragmentation at the Fermilab Tevatron. Fig. 1(a) is one of the 
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Feynman diagrams of the gluon-gluon fusion processes for the DP interactions, and Fig. 1(b) 
for the SP interactions. In our calculations of the DP cross section for the large px ip'ip 
production, we only consider the relevant single J/ip production at large px via color-octet 
gluon fragmentation. For the double J/-0 production from the SP interactions, we also 
calculate the contributions only from color-octet gluon fragmentation [17], which dominate 
over the color-singlet contributions both from the leading-order processes [ID] and other 
fragmentation processes |12|] at large px due to the same reason as pointed out in fl4[ . 

In the DP interaction processes, the two partonic interactions occur independently of 
. Toll . So, the cross section for 'ifj'ijj production from the DP processes can be 


each other 


related to the cross section for single J/-0 production from the SP processes by 

(?spW 


CTDp(lp'tp) 


2a, 


eff 


( 3 ) 


where the effective cross section a e ff represents the possible correlation effects of the parton 
distributions in the proton (antiproton). If parton correlations are negligible, 2a e // should 
approximately be equal to the total inelastic cross section of 44m6 at the Tevatron []li|. This 
implies a e // ~ 22 mb. Parton correlations tend to reduce the effective cross section (i.e., 
increase the double parton scattering cross section) relative to the uncorrelated case. In the 
following calculations, we use the value a e ff = 14.5mb according to the measurement by the 
CDF ||. In the literature, some modifications to the above formula have been introduced for 
the calculation of the DP cross section, which may more correctly represent the correlations 
between the two partons in one proton such as the energy-momentum conservation effects 
[7j.p~0f] . In our calculations, we do not consider these modifications, because they only cause 
a little change to the total cross section and can then be neglected [7[|1CS . 

To calculate the single J/iJj production rate, we consider qq/gg —> gg and q(q)g —■> q(q)g 
subprocesses, and then gluon fragmentation to J/'ijj. In gluon fragmentation, the input 
parameters are taken to be 


m c = 1.5 GeV, a s (2m c ) = 0.26, (O^ ( 3 Si)) = 0.0106GeD 3 . 


( 4 ) 


The value of the color-octet matrix element ( 3 Si)) follows the fitted value in p0| by 

comparing the theoretical prediction to the experimental data at the Tevatron. We use the 
MRS (A) parton distribution functions [j2l] to generate the production cross section, and set 


4 




































the renomalization scale and the factorization scale both equal to the transverse momentum 
of the fragmenting gluon fi = pr(g) ~ Pt{ VO- A pseudorapidity cut of \r]{ij))\ < 0.6 was 
also performed on the produced ips. We obtain the integrated cross section for single J/ijj 
production (over pr) as a function of the minimum pr(V0- In the calculations of gluon 
fragmentation, we also include the contributions from the \ c and ip ' feeddowns through 
g Xc and g —> ip' followed by y c —> -07 and ip' —> ipX. This means that the calculated 
cross section is for the prompt production (excluding the contributions from b decays). The 
feeddown contributions give the same pt distribution of J/ip and increase the total rate by 
a factor « 1.6. We estimate this factor from the measured fraction of direct production in 
the prompt production [|l3j (which is 64%). The leptonic decay branching ratio Br(J/'i/j —> 
p + /i~) = 0.0597 is also multiplied in the cross section. 

Substituting the integrated cross section for single J/ip production into Eq.([|), we can 
estimate the cross section for double J/ip production coming from the DP processes. The 
final results are plotted in Fig.2 and Fig.3 for the experiments at the Fcrmilab Tevatron and 
at the CERN LHC respectively. The curves represent the integrated cross sections for double 
J/ip production (over p T ) a(pp(p) —> ipip + X ) x Br(J/ip —> /i + /i - ) 2 as the functions of the 
minimum p?( -0). For comparison, we also plot the contributions from the SP interactions 
via double gluon fragmentation within a single parton-parton scattering. We calculate these 
contributions by using qq/gg —> gg subprocesses followed by the two gluons fragmentation 
into two J/ij )s 0- Considering the dominance of the gluon-gluon fusion subprocess in single 
J/ip production, as a rough estimate, we can also write the SP ipip production cross section 
as flTjJ 

1 


^SP(#) ~ 7WSP(V0 X fg^, 


(5) 


where f g is the gluon fragmentation probability to J/ip at large Pt- In these two figures, 
the solid lines represent the contributions from the DP interactions, and the dotted lines 
correspond to the contributions from the SP interactions. 

As for the DP processes, suitable kinematic cuts can be used to detect the signals. The SP 
processes, gg —> ijnp produce the J/ip pair back-to-back in transverse momentum, whereas 
the DP processes produce the J/ip pair unrelated to each other but balanced in p t by a hard 
gluon (quark) jet. 

From Fig.2, we can see that at the Tevatron the DP cross section is smaller than the SP 
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cross section in all pr(min) region. Prompt single J/0 production rate at the Tevatron has 
been measured with pr(0) > 4 GeV and |t/( 0)| < 0.6, and the integrated cross section is 

a(pp —■> J/0 + A") x Br(J/ip —> /i + /i~) ~ 24nb. (6) 


Substituting the above value of the single J/0 production cross section into Eq.(|3|), we obtain 
the double J/0 production cross section from the Double Parton processes, 

a(pp -> 00 + X ) x Hr(0 -»/i + /Q 2 = 20fb = 6 

cr (pp —> 0 + X) x Br(ip —> p + p~) 24nb 

which shows that there will be about one 00 event among every 10 6 single 0 events. The 
fraction of double J/0 events from the DP processes to single J/0 events is an order of 
magnitude smaller than that from the SP process (T^j (where it is 7.6 x 1CP 6 ). The cross 
section for double J/0 production from the DP processes is estimated to be about 20fb with 
Pt(0) > 4 GeV and |7/(0)| < 0.6. This indicates that the DP -0-0 signal is at the edge of the 
detectability of the Tevatron at present (considering the integrated luminosity of 100pb _1 
now accumulated by each of the Tevatron detectors and the possible inclusion of both the 
J/0 —> / u + /i“and J/0 —> e + e _ modes), and future increases in luminosity could possibly 
make it measurable. 

However, the DP cross section for double J/0 production is proportional to the square of 
the single J/0 production cross section (see Eq.(|)), whereas the SP cross section for double 


J/0 production is proportional to single J/0 production cross section (see ||T7j). The relative 
importance of the DP contributions to 00 production against the SP contributions will be 
changed with the single J/0 production rate. At high enough energies, the DP contributions 
may dominate over the SP contributions because the single J/0 production rate increases 
with the collider energy. To see this more clearly, we give the explicit expression for the ratio 

crsp(00) (Teff X f g ^ 

which is obtained from Eqs.([|) and (|5]). At the Tevatron the ratio is about 0.11. But 
the ratio will increase as erspi' 1 !’) increases with the collider energies. At the LHC, we find 
the single J/0 production rate will be over an order of magnitude higher than that at the 
Tevatron (see also |22|). So, the DP contributions to the double J/0 production will be more 
important at the LHC, which is shown in Fig.3. This figure shows that in the lower pr(min) 
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region (i.e., pp(min) < 7 GeV) the DP contributions dominate over the SP contributions. 
For pp(min) = 5 GeV, the DP cross section is 

(r D p{i>i’) = 1.45pb, (9) 


while the SP cross section is 


asp^ip) — 0.63pb. 


( 10 ) 


In the above calculations, we choose the effective cross section at the LHC as the same as 
that at the Tevatron. In fact, the value of a e ff may change its value while the collider 
energy y/s increases. However, we notice that the total cross section increases slowly as y/s 
increases, (e.g., a tot ~ lOOmb at the LHC while a tot ~ 80mb at the Tevatron 0). So, we 
expect that cr e ff will not change much at the LHC. 

For a typical integrated luminosity ~ 10 4 pb _1 at the LHC, we would expect the order 
of 10 4 events of DP , 0'0. This indicates that we can detect the DP signals and also measure 
the cr e ff to investigate the possible correlations between partons in the proton. 

It should be emphasized that the color-octet production mechanism is of crucial impor¬ 
tance to the double J/ij} production from the DP interactions at the Fermilab Tevatron and 
the CERN LHC, just as it is to the single J/^ production at Tevatron M, 


Within the 


color-singlet model, the DP contribution is much smaller than that from the SP interaction, 
i.e., <T£ip('0'0) •C u S p('0'0) ©]. However, after including the color-octet mechanism, the 
DP contribution can dominate over the SP contribution for lower pp(min) region at the 
LHC. So, double J/if; production at large pp can also provide another important test for the 
color-octet production mechanism. 

In conclusion, we have calculated double J/%b production at hadron colliders. We find 
that the DP production is at the edge of the detectability of the Tevatron at present, 
and at the LHC the DP cross section will dominate over the SP cross section in the lower 
Pr(min) region (i.e., pp(min) < 7 GeV). We also find the new production mechanism, i.e., 
the color-octet mechanism is of crucial importance to the double J/-0 production at high 
energy hadron colliders. Therefore, the measurement of the double J/-0 production rate 
would provide an important test for both the DP scattering and the color-octet production 
mechanism. 
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Figure Captions 


Fig.l. The typical Feynman Diagrams for double J/i^> production at large pr via color-octet 
gluon fragmentation at the Fcrmilab Tevatron. One of the diagrams of the gluon-gluon 
fusion subprocesses for (a) the double parton interactions, and for (b) the single parton 
interactions. 

Fig.2. The integrated cross section of t/u/ ; production x Brfy —> /j, + p~) 2 for p T > p T (min) 
as a function of minimum pri^) at the Fcrmilab Tevatron. The dotted line represents the 
contributions from the SP processes via double gluon fragmentation, and the solid line 
corresponds to contributions from the DP processes. 

Fig.3. The integrated cross section of ipip production at the CERN LHC. Here the curves 
are defined as those in Fig.2. 
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